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ABSTRACT 


Oil  exploration  in  water  deeper  than  1000  feet  has  motivated  the  development 
of  drilling  platforms  with  dynamic  responses  superior  to  those  of  the  conven¬ 
tional  jacket  structures.  One  of  the  most  promising  of  these  is  the  tension.,  , 
leg  platform  (TIP),  a  semi -submersible  platform  held  in  place  by  tension  members 
connected  to  the  ocean  bottom. 

This  thesis  evaluates  the  wave  force  transfer  function  for  the  TLP  in 
heave,  pitch,  and  roll  motions.  This  is  accomplished  by  assuming  the  transfer 
function  for  a  single  leg  and  then  extending  that  result  to  include  the  effect 
of  phase  differences  between  forces  on  multiple  legs,  including  the  effect  of 
wave  spreading. 

Cosine  squared  directional  wave  spreading  is  discussed  and  then  dismissed 
in  favor  of  a  Gaussian  distribution  of  wave  direction.  Using  the  Gaussian 
distribution  it  is  shown  that  when  the  wave  spreading  has  a  standard  deviation 
greater  than  0.5  radian  the  exciting  force  transfer  functions  in  heave,  pitch, 
and  roll  converge  rapidly  to  that  computed  for  seas  uniformly  distributed  over 
all  incidence  angles,  y 
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ABSTRACT 

Oil  exploration  in  water  deeper  than  1000  feet  has 
motivated  the  development  of  drilling  platforms  with 
dynamic  responses  superior  to  those  of  the  conventional 
jacket  structures.  One  of  the  most  promising  of  these 
is  the  tension  leg  platform  (TLP) ,  a  semi-submersible 
platform  held  in  place  by  tension  members  connected  to 
the  ocean  bottom. 

This  thesis  evaluates  the  wave  force  transfer  function 
for  the  TLP  in  heave,  pitch,  and  roll  motions.  This  is 
accomplished  by  assuming  the  transfer  function  for  a 
single  leg  and  then  extending  that  result  to  include  the 
effect  of  phase  differences  between  forces  on  multiple 
legs,  including  the  effect  of  wave  spreading. 
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Cosine  squared  directional  wave  spreading  is  discussed 
and  then  dismissed  in  favor  of  a  Gaussian  distribution  of 
wave  direction.  Using  the  Gaussian  distribution  it  is 
shown  that  when  the  wave  spreading  has  a  standard  deviation 
greater  than  0.5  radian  the  exciting  force  transfer 
functions  in  heave,  pitch,  and  roll  converge  rapidly  to 
that  computed  for  seas  uniformly  distributed  over  all 
incidence  angles. 
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I .  Introduction 


Present  day  platforms,  in  water  depths  from  10  to  1000 
feet,  are  primarily  of  welded  steel  jacket  construction.  As 
the  present  limit  of  about  1000  feet  was  approached  many 
problems  in  fabrication  and  installation  were  encountered, 
indicating  the  need  for  a  shift  to  a  new  technology.  More 
significantly,  with  inceasing  water  depth  the  natural  per¬ 
iods  of  jacket  platforms  increase  from  about  one  second  in 
100  feet  of  water  to  approximately  four  seconds  in  1000  feet 
of  water.  There  is  significant  wave  energy  with  a  period  of 
four  seconds  and  dynamic  amplification  becomes  a  problem. 

The  expected  life  of  the  structure  is  shortened  due  to  high 
cycle  -  low  stress  fatigue  damage. 

One  alternative  structure  is  a  tension  leg  platform 
(TLP) .  A  TLP  (see  Figure  1-1)  is  a  bouyant  structure  held 
in  position  by  tension  members  secured  to  the  ocean  bottom. 
Because  the  platform  is  not  rigidly  connected  to  the  bottom 
its  natural  modes  of  vibration  are  significantly  different 
from  a  jacket  structure  of  similar  height.  A  TLP  in  a  water 
depth  on  the  order  of  1000  feet  will  have  natural  periods 
in  heave,  pitch,  and  roll  of  about  two  seconds,  well  out  of 
the  region  of  high  wave  energies. 
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Since  the  TLP  has  only  recently  been  accepted  as  a 
potential  alternative  for  deep  water  drilling  there  are  few 
studies  of  its  dynamic  response  to  wave  forces.  Such  preli¬ 
minary  studies  are  necessary  to  provide  information  for  the 
design  of  the  operating  platforms.  Response  prediction  re¬ 
quires  specification  of  the  modal  exciting  forces,  for  example 
in  pitch,  roll,  and  heave.  For  low  everyday  sea  states  the 
exciting  forces  and  structure  responses  can  be  considered  to 
behave  in  a  linear  fashion.  Linearity  of  excitation  and  re¬ 
sponse  suggests  that  a  frequency  domain  solution  technique 
is  appropriate;  especially  when  one  considers  that  ocean 
waves  are  a  random  process  most  often  described  by  a  wave 
amplitude  spectrum.  This  thesis  presents  an  estimate  of  the 
modal  exciting  force  spectra  in  heave,  pitch,  and  roll  in  a 
way  which  accounts  for  structure  geometry  and  wave  spreading. 
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II.  Problem  Statement 


The  modal  exciting  forces  on  a  tension  leg  platform  (TLP) 
will  be  dependent  upon  the  frequency  and  the  angles  of  inci¬ 
dence  of  the  incoming  waves.  For  some  angles  of  incidence 
(at  a  given  frequency)  the  forces  on  the  various  legs  will 
act  together  to  produce  a  larger  total  force  or  moment.  For 
other  angles  of  incidence  the  phase  differences  will  cause 
the  forces  or  moments  exerted  on  the  various  legs  to  cancel, 
tending  to  produce  a  smaller  total  force. 

In  this  report  the  total  modal  force  of  the  platform  in 
heave,  pitch,  and  roll  will  be  divided  into  two  parts.  The 
first  will  be  the  force  exerted  upon  a  single,  isolated  leg 
by  waves  of  the  given  frequency.  The  second  will  be  a  geo¬ 
metric  term  to  take  into  account  the  effect  of  the  phase 
differences  between  forces  on  each  of  the  legs.  This  second 
term  will  be  specifically  calculated  for  a  square  TLP  noting 
that  the  method  could  be  used  for  any  shape. 

The  calculations  so  far  imply  a  unidirectional  wave 
amplitude  spectrum.  Rarely  is  this  observed  in  practice. 

A  reasonable  wave  spectrum  can  be  approximated  as  being 
centered  about  one  bearing  from  the  structure  and 
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diminishing  to  zero  at  bearings  90  either  side  of  the 
central  bearing.  Two  types  of  wave  spreading  functions 
will  be  discussed.  The  first  is  the  cosine  squared  spread¬ 
ing  law  which,  although  widely  used,  allows  little  if  any 
leeway  in  varying  the  spreading  to  fit  actual  measurements. 
The  second  is  the  Gaussian  distribution  which  does  allow 
continuous  variation  of  spreading  simply  by  changing  the 
standard  deviation. 

The  modal  wave  force  spectra  upon  the  structure  will 
be  the  integral  over  all  incidence  angles  of  the  directional 
wave  amplitude  spectra,  weighted  appropriately  by  the  plat¬ 
form  heave  force,  roll  or  pitch  moment  transfer  functions. 
The  resulting  spectra  will  be  normalized  in  this  thesis  by 
the  results  obtained  when  a  uniformly  spread  wave  spectrum 
is  used.  This  is  done  to  demonstrate  the  effect  of  wave 
spreading  on  the  heave,  pitch,  and  roll  exciting  force 
spectra . 
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III.  Formulation  of  Equations 

A.  Platform  transfer  function 

1.  Single  leg  transfer  function.  The  single  leg  wave 
amplitude  to  vertical  force  transfer  function#  rs(&>)»  Is  a 
function  only  of  wave  frequency.  Because  the  leg  has  circu¬ 
lar  cross  sections  there  is  no  dependence  on  $#  the  angle  of 

incidence.  This  thesis  assumes  that  rs(ti»)  is  available# 
either  experimentally  or  through  another  approximation  method. 

2.  Square  platform  transfer  function.  Consider  plane 
progressive  waves  of  frequency  u  coming  from  an  angle  g  off  the 
platform  axis.  See  Figure  XXX- 1.  Bach  leg  is  thus  subject  to 
waves  of  the  same  frequency  but  different  phase.  For  deep 
water  waves  the  wavenumber  is  k  ■  h  ■  u2 .  Thus  the  wave- 

X  g 

length  is 


The  phase  shift  in  one  wavelength  is  2ir.  The  phase 
shift#  4^#  due  to  a  distance#  1^,  along  the  axis  of  prop¬ 
agation  is  thus 


Consider  the  center  of  the  platform  as  the  point  of  sero 
phase  shift.  Let  d  be  the  leg  spacing  and  a  -  d//2  be  the 
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Figure  III-l  Regular  Waves  Incident  on  TLP 


distance  from  the  center  of  the  platform  to  the  center  of  any 


leg. 


The  distance  causing  a  phase  shift  for  leg  1  is 


1^  »  a  cos(-  -  0) .  Por  leg  2  it  is  lj  ■  -a  sin(j  -  0) . 

Similarly  for  legs  3  and  4  the  lengths  are  13  -  a  sin( 

and  1.  -  -a  cos (7  -  0).  From  Equation  XXX-1  the  phase 
4  4 


shifts  are: 


4  m  2  a  cos (7  -  0) 

i  g  4 


a  ■»«  a  sin (7  ~  0) 
a  g  4 

4  -  S*a  sin (7  -  0) 
a  g  4 

4  a  cos (7  -  0) 
%  g  4 


III-2 


Mote  however  that: 


cos(*  -  0)  -  cos (2) cos (3)  ♦  sin(|)sin(0) 

■  (/2/2) {cos (0)  +  sin(0) } 

sin(J  -  0)  ■  sin (J) cos (0)  -  cos (J) sin (8) 

■  (/2/2) (cos (0)  -  sin(0)> 


Letting  T  ■  «*d/g  and  recalling  a  -  (/2/2)d  then 
Equations  XXX-2  become: 
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t  |  (cos(8)  +  sin(8)> 

1  z 

*  — S  {cos (8)  -  sin(8) ) 

2  2 

•  I  (cos(6)  -  sin (8) ) 

+  (costs)  +  sin(8) )  m-3 

3.  Heave.  The  transfer  function  for  heave  for  the 
total  structure,  rht,  is  the  sum  of  the  contributions  of 
each  leg.  Taking  phase  into  consideration  each  leg  contri¬ 
butes  rs  exp(i(<ut+$i)  }.  Thus  the  quantity  of  interest,  the 
square  of  the  magnitude  of  is 

lrhtl2  -  |ra|  *  |eA*l+  el*a+  el*»+  *L*'\*. 

Incorporating  the  identity  e*^  ■  cos($)  +  i  sin($) , 
summing  real  and  imaginary  parts,  and  taking  the  magnitude 
squared  of  the  exponentials  this  becomes 

|rhtl*  ■  |r#|*  ({cos($j)+cos($a)+cos($j)+cos($%) )* 

♦{sin($i)+sin($2)+sin(4,)+8in(4%) }a)  ZZZ-4 

Mote  however  from  Equation  ZZZ-3  that  4*  -  -4t  and 


♦a  ■ 


Thus  we  have 


cos(4x)  »  cos(4„)  sin(4x)  ■  -sin(4*) 

cos(42)  »  cos(4x)  sin(42)  ■  -ain(42) 

and  Equation  ZZZ-4  becomes 

|rht|*  *  |r,j*  (2  cos($x )  +  2  cos(4t)}a 

|rht|*  -  4  |rja  (coa(4x)  ♦  cos(4j)}a.  ZZZ-5 

By  replacing  4X  and  4S  with  their  valuea  from  Equation 
ZZZ-3  and  netting  ■  J rhtl  a/|rs|  *  becomes 

Bjj  *  4  (coa{|(cos  0  +ain  0  ) }  ♦  cos {|( cos  0  -sin  0  )})a 

f  y  T  - 

Rjj  *  4  (cost-  cos  0)coa(^  sin  0)  -  sin(j  cos  0)sin(^  sin  0) 

t  coa(j  cos  0)cos(^  sin  0)  +  sintj  cos  B) sin (|  sin  0)}1 

R»*  4  {2  cos(|  cos  0) cos ( “  sin  0))a  ZZZ-6 

n  2  2 

A  form  often  more  suitable  for  computer  work  is  obtained 
by  applying  the  identity  cos*  A  *  ij  +  %  cos  2A  yielding 

|rht|a  »  4  (1+  coa(T  cos  0))  (!♦  coa(T  sin  0))  | | * 

ZZZ-7 

which  can  be  implemented  by  recalling  that  T*  «ad/g. 


16 


4.  Horizontal  force  assumptions.  The  roll  and  pitch 
moments  experienced  by  a  TLP  will  be  dependent  upon  horizon¬ 
tal  as  well  as  vertical  forces.  In  this  report  it  has  been 

assumed  that  the  total  moment  producing  effect  of  the  hori- 

% 

zontal  force  is  small  compared  to  the  vertical  forces.  This 
is  perhaps  not  always  the  case  but  is  assumed  in  this  discus¬ 
sion  because,  while  the  horizontal  wave  forces  may  be  large, 
they  are  distributed  along  the  legs,  both  above  and  below  the 
center  of  rotation  of  the  platform.  The  force  on  the  leg 
near  the  surface  is  large  but  with  a  small  moment  arm.  The 
forces  below  the  center  of  rotation  will  be  smaller  but  will 
extend  further  down  the  leg  and  thus  have  a  large  moment  arm, 
producing  a  moment  similar  in  magnitude  but  opposite  in  sign, 
and  will  tend  to  cancel  the  other.  Thus  the  moment  producing 
effects  of  the  horizontal  forces  are  assumed  to  be  signifi¬ 
cantly  smaller  than  those  of  the  vertical  forces.  When  the 
location  of  the  center  of  rotation  and  the  roll  and  pitch 
natural  periods  are  known,  as  they  would  be  for  any  particu¬ 
lar  TLP  design,  the  horizontal  contribution  to  the  roll  and 
pitch  moment  could  be  calculated.  It  can  be  shown  that  the 
phase  angles  encountered  when  adding  up  the  moments  generated 
by  the  horizontal  forces  on  each  of  the  four  legs  would  be 
the  same  as  calculated  previously  for  heave.  Therefore  the 


17 


results  computed  for  heave  will  serve  as  a  check  on  the  qual- 
itative  behavior  of  the  role  of  horizontal  forces  in  the 
generation  of  roll  and  pitch  exciting  forces. 

5.  Pitch  and  roll.  Pitch  and  roll  for  a  square  plat¬ 
form  are  exactly  symmetric.  Thus  only  roll  will  be  discussed 
here  with  the  understanding  that  pitch  could  be  obtained  by 
shifting  the  coordinates  by  90  degrees.  Taking  the  direction 
of  positive  roll  moment  to  be  as  shown  in  Figure  IIX-1  the 
total  transfer  function  for  roll,  squared,  would  be 

lrrtla  *  |T  l*  | (d/2) e1* 1 ♦ (d/2) e1*  *♦ ( -d/2) e1*  * ♦ ( -d/2) e1* * | * 

*  lr»l*  f  | ei^l+ei^2-ei**-ei^'*  |  2 

Let  -  |rrt| */|rs| a.  Incorporating  the  identity 
e*4  ■  cos  ♦  i  sin  $  ,  summing  the  real  and  imaginary 

parts,  and  taking  the  magnitude  squared  this  becomes 

£2 

Rr  *  |  {(cos  cos  ♦j-  cos  ♦j-  cos 

♦(sin  ♦,♦  sin  sin  sin  ♦„)*)  III-8 

However,  from  equations  III -3  4*  ■  -♦ j  and  +  s  ■  -$2  . 
Thus  cos  ■  cos  ,  cos  ♦  ,  -  cos  $2  #  sin  ■  -sin  , 
and  sin  -sin  $2  .  Therefore  the  cosine  term  in 
Equation  III-8  reduces  to  zero  and  the  sine  term  becomes 
(2  sin  ♦  2  sin  +a)**  leaving 
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d2  (  sin  $|+  sin  $2) 2 . 


Incorporating  the  values  from  Equations  III-3  this  becomes 

Rr  -  da  (sin{|(cos  8  +sin  6  ) }+sin{-|(cos  8  -sin  6  )))* 

R  *  dl  {sin(|cos  B)cos(|sin  3)  +  cos(Icos  8 ) sin (f sin  0) 
t  2  2  2  2 

+  sin(-|cos  8)cos(|sin  0)  +  cos(-|cos  8)sin(^sin  8)) 2 

R  »  d2  {sin(-cos  8)cos(-;sin  0)  +  cos (f cos  0)sin(£sin  0} 
r  2  2  2  2 

-  sin(|cos  0)cos(|sin  0)  +  cos(^cos  8)sin(|sin  0)>2 

R  ■  4  da  cosa(|cos  0) sin2 (fsin  8)  III-9 

r  2  2 

By  applying  the  identities  cos2A  »  h  +  **  cos  2A  and 
sinaA  m  h  -  h  cos  2A  this  yields  a  form  more  suitable  for 
computer  work. 

|  Trt |  2  *  da  {1  ♦  cos(T  cos  8)  Hi  -  cos  (T  sin  8)}  jrs|2 

III-10 

As  mentioned  previously  pitch  is  singly  rotated  by  90 
degrees.  Thus  cos  8  becomes  sin  $  and  sin  0  becomes  -cos  0  . 
The  transfer  function  for  pitch  for  a  square  platform  is  thus 

Jr  . J*  ■  da  {1  -  cos(T  cos  0)}{1  ♦  cos(T  sin  0))  |r_|*. 

pt  » 

III-ll 
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B.  Wave  spectrum. 

The  wave  spectrum  incident  on  the  platform  varies  with 
both  frequency  and  direction  and  is  denoted  Sn(oj,B,8#). 

0,  is  the  center  of  an  arc  of  180°  from  which  the  waves 
are  coming. 

For  simplicity  of  calculations  it  is  often  assumed,  as 
it  shall  be  in  this  report,  that  the  total  wave  spectrum 
can  undergo  separation  of  variables.  Thus  we  can  express 
the  spectrum  as 

S  (a»,8,B.)  -  S  (at)  F(0,B  ) 

n  •  n  * 

where  the  entire  frequency  dependence  is  contained  in 
Sn(o>)  and  the  entire  angular  dependence  is  in  F(0,0a). 

This  section  deals  with  approximations  for  the  angular 
dependence ,  F ( 8 , 6  0 ) . 

1.  Cosine  squared  wave  spectrum.  An  often  used  method 
for  describing  the  angular  spreading  of  waves  is  the  cosine 
squared  wave  spectrum  (see  Figure  III-2) .  It  has  the 
advantage  that  it  vanishes  for  values  of  8  equal  to  plus  and 
minus  ir/2.  The  formula  for  the  cosine  squared  spectrum  is 

Sn(a>,e,B,)  -  C  p(8,0t)  Sn(a>)  *  -tr/2  <  8-8,  <  v/2 


where  8,  is  the  center  of  the  wave  spectrum  and  C  is  a  constant 


used  to  normalize  to  unity  the  area  between  -ir/ 2  and  it/ 2. 
It  is  derived  as  follows: 


ir/2 

1  *  /  C  cos2 (8)  d8 
-ir/2 

1  */2  * 
i  -  2  /  cos2  (6)  d0  «  2  <5> 

C  4 


Thus  the  equation  for  the  cosine  squared  spectrum  is 
Sn(e,0,0a)  -  2  Sn(«)  cos2 (0-0,)  ;  <  (0-0,)  <  | 

111-12 

This  spectrum  suffers  from  one  major  disadvantage. 

It  cannot  allow  for  variations  in  the  directionality  of  the 
waves.  This  problem  can  be  eased  somewhat  by  considering 
cosine  to  the  fourth,  to  the  sixth,  and  so  forth  but  this 
results  in  only  step  changes  in  the  directionality.  Methods 
have  been  designed  to  provide  continuous  changes  but  they  are 
somewhat  cumbersome  and  comparison  is  not  apparent. 

2.  Gaussian  wave  spectrum.  The  method  used  in  this  re¬ 
port  to  achieve  continuous  changes  in  directionality  is  to 
consider  the  waves  to  be  spread  about  the  central  angle  0(  in 
a  Gaussian  fashion  with  standard  deviation  o.  The  reader  not 
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familiar  with  Gaussian  (or  normal)  distributions  may  wish  to 
refer  to  any  good  probability  text  for  an  explanation  more 
detailed  than  the  following. 

The  probability  distribution  of  a  Gaussian  process  with 
mean  8«  and  standard  deviation  o  is  described  by 

p(8,ft0)  ■  — - —  exp(-(B-8#) a/(2o2) } 

/2ir  a 

A  spectrum  of  waves  defined  by  a  Gaussian  process  would 
thus  be 

SnU,0,Bo)  *  C  p(8,8,)  Sn(«)  ;  -ir/2<8-8,<ir/2  III-13 


where  again  C  is  a  constant  to  normalize  to  unity  the  area 
between  -ir/2  and  ir/2.  It  is  derived  as  follows,  referencing 
8t  to  zero: 

ir/2 

1  -  /  C  p(8,6«)  dS 

-ir/2 

ir/2 

-  ■  /  p(S,8» )  d8 

c  -tt/2 


This  integral  is  evaluated  using  the  properties  of  the 
Gaussian  distribution  and  yields 


C 


)-l 


x 

/  exp(-ta/2)  dt 


where  t(x) 
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is  known  as  the  normal  distribution  function  and  is  tabulated 
in  probability  texts  and  math  tables.  Thus  C  is  a  function 
of  the  standard  deviation  and  shall  hereafter  be  denoted  C(a). 

Unlike  the  cosine  squared  spreading  the  Gaussian  distri¬ 
bution  does  not  vanish  for  0  equal  to  plus  and  minus  ir/2. 

This  is  why  C(o)  varies  with  o.  It  must  account  for  the 
cropping  of  the  "tails"  of  the  distribution  for  1 8 1  >  ir/2. 

Some  simplifying  approximations  can  be  made  however. 

As  can  be  seen  from  Figure  III- 3  when  a  reaches  a  value  of 
2  or  larger  the  angular  spectrum  approaches  a  constant 
value.  That  value  is  1/ir  (in  order  to  achieve  unity  area) . 
Thus  Equation  III-13  becomes 

S_(u,0,0,)  *  3  S  (w)  *  ””  <  B-0o  K  r  for  a>2. 

14  n  n  2  2 

For  o  <  1  the  area  in  the  tails  diminishes  rapidly 

leaving  almost  all  of  the  area  between  -ir/2  and  ir/2.  As 

can  be  seen  in  Figure  III-4  C(o)  rapidly  approaches  unity. 

For  0*1.0  01.13.  For  0*0.9  01.09.  So  for  a  <  1 

Equation  III-13  can  be  approximated  by 

SnU,0,0o)  *  S  (<a)  {-JL — }  exp{~(^~f^-}  . 

n  n  /2ir  o  20 

For  values  of  0  between  1  and  2  the  entire  exact  formula 


must  be  used 


•to..,.,)  -  ..to)  .xp(^f^)  m-l« 

n  20 

3.  Comparison  of  cosine  squared  with  Gaussian.  Since 
many  people  are  familiar  with  cosine  squared  spreading  it 
is  useful  to  discuss  for  what  values  of  a  the  Gaussian  would 
approximate  the  cosine  squared.  Although  the  two  functions 
are  quite  distinct  there  are  several  possible  methods  to 
compare  them. 

The  normalized  cosine  squared  distribution  has  a  value 
of  2/ir  *  0.637  for  8  ■  B».  For  Gaussian  spreading  this  would 
be  C(a)  /  (/2v  a)  at  8  *  8« .  If  a  is  assumed  less  than  1 
then  C(o)  3  1.  Equating  the  two  values  yeilds  a  *  0.627. 

When  B  «  8,  +  ir/4  the  cosine  squared  distribution  is 
half  of  its  maximum  value.  This  situation  can  be  matched 
with  a  Gaussian  yielding  still  another  "equivalent"  a  for 
the  cosine  squared. 

C (at)  .-(tt/4)  * .  CU) 

«*pt  -  ii«— >  -  <0-51  71T7 

expC^r™ )  "  (0.5) 

2cz 


a  m 


f-U/4)» 

1  ln(0. 5)  1 


a  *  0.667 
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The  two  values  obtained  are  approximately  equal.  Thus  the 
Gaussian  distribution  with  c*0.65  is  very  similar  to  the  more 
familiar  cosine  squared  distribution.  A  smaller  cr  means  a  more 
unidirectional  wave  front  and  a  larger  o  means  more  spreading. 

C.  Platform  exciting  forces. 

Each  differential  element  of  the  spectrum  Sn(a,8,Sa)  dB 
induces  mean  square  force  or  moment  upon  the  platform  through 
the  appropriate  platform  modal  force  transfer  function, 
magnitude  squared,  |rt(o,8)|*.  Thus  the  total  force  or 
moment  spectrum  is  given  by 


H(«d,8« )  -  /  |T  U,8)|a  S  (0,8,8,)  dB  III-15 
-w  t  n 

It  is  important  to  note  that  this  modal  force  spectrum 
is  dependent  upon  8, ,  the  center  of  the  directional  spectrum. 

1.  Square  platform,  heave  force.  The  heave  force 
spectrum  for  a  square  platform  in  normally  distributed 
waves  is  derived  by  substituting  Equations  III-7  and  HI-14 

into  Equation  HI-15  yielding: 

B,+*/2 

H K(o, 8,)  ■  /  |T  |*  4{l+cos(T  cos  B)Hl+cos(T  sin  8))* 

B,-ir/2  * 


dB 
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I 


Hh(«a#8,) 


B(^/2  »a  a  i  1 

/  {l+cos(T  cos  8)  Hl+cos  (T  sin  6> )  exp{  )  dp  III-16 

8, -s/2  2  0 


Sotting  Xh(u,P  )  equal  to  the  integral  in  Equation 
III-16  and  expanding  the  integrand  yields 


Pi-fs/2 

I  (a»,P,)  *  /  (1  +  cos(T  cos  8)  ♦  cos (T  sin  P)  + 

h  8, -s/2 

cos(T  cos  P)cos(T  sin  8))  exp{  ^a**'  ^  III-17 


where  Hh<«,8,)  -  4  IrJ*  Sn(«)  IhU,S,)  • 

A  dosed  form  solution  of  Equation  III-17  does  not 
exist.  Thus  a  numerical  integration  was  necessary 
(  see  Section  IV  ) . 

2.  Square  platform,  roll  moment.  Similarly  the  roll 
moment  spectrum  for  a  square  platform  in  normally  distributed 


waves  is  derived  by  substituting  Equations  III-10  and  XXJ-14 

,/ 

into  Equation  IXX-15  yielding: 


B,+»/2 

H  (<o,3»)  *  /  d*  (l+cos(T  cos  8)Hl-cos(T  sin  8)}  |T  |  ** 
S.-T/2  * 

S  (m)  {-—  i  ■  1  (jmL — }  exptrJ8.-6.*)*)  d6 

n  2  ♦  (*  )-l  /2n  a  2 

la 


H  (ai#8,)  -  d*  |r  |a  S  («)  x 

r  s'  n  / 2ir  o 


— (B— g.\* 

/  {l+cos(T  cos  8)){l-cos(T  sin  8)}  exp{ — — )  df}  III-18 
8,-ir/2 


Setting  I  (<u,Sa)  equal  to  the  integral  in  Equation  111-18 

and  expanding  the  integrand  yields 
Bt+ir/2 

I  U,80)  »  /  {1+cos (T  cos  8) -cos (T  sin  8)- 

8,-ir/2 


cos(T  cos  8)cos(T  sin  8) >  exp{ 


-(8~Bb ) * 

2o* 


}  d8 


III-19 


where  Hr(<u#B,) 


d*  |r,|* 


sn«“) 


C(q) 
/ 2n  a 


Ir(o»»8»)  . 


Again  a  closed  fora  solution  does  not  exist  and  the  same 
numerical  integration  technique  was  applied  (  see  Section  IV) . 
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3.  Square  platform,  pitch  moment.  As  mentioned 
previously,  for  a  square  TLP  the  pitch  moment  is  simply 
rotated  90  degrees  from  the  roll  moment.  Thus 


-  a*  |r,|*  snu>  tw  i  !.,(,) 


where  St+ff/2) 


4.  Normalisation .  In  order  to  compare  the  effect  of 
the  wave  spreading  for  different  configurations  a  method 
of  normalisation  was  developed.  The  normalisation  factor 
was  taken  to  be  the  force  or  moment  frequency  spectrum 
resulting  from  waves  of  equal  amplitude  from  all  directions 
in  one  half  of  the  plane.  Thus  SnU,8,Ba)  equals  (1/v)  s„u> 
for  (B«-ir/2)  <  0  <  (0«<Mr/2)  and  sero  elsewhere.  The 
normalisation  frequency  spectrum  becomes 

0o+ir/2 

H  U, 8,)  -  /  |r  U,B)|*  (1/ff)  S  U>  d8  . 

0  S,-ir/2 

Because  of  the  symmetry  of  a  square  platform  and 
the  uniformity  of  this  spreading  it  makes  no  difference 
where  the  center  of  the  spectrum  is.  For  convenience 
0,  was  taken  to  be  sero  and  thus  the  normalisation  spectrum 


becomes 


1  »/2 

■  Z  sn(w)  /  |rt(ir,8)  |a  d6 

*  -ir/2  c 


IXI-20 
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which  is  simply  the  mean  square  value  of  r*  times  S  (to)  . 

c  n 

For  heave  and  roll  they  are  as  follows 


4  ir/2 

(o») »—  S_  (ci>)  |r|*  /  (l+cos(T  cos  B )  }  ( 1+cos  (T  sin  B) }  dB 

0h  *  "  •  -ir/2 

III-21 


A*  */2 

Ho -(*)■"  S  {«)  Jr  I1  /  { l+cos (T  cos  B)){l-cos(T  sin  B) )  dg 
or  r  n  s  _t/2 

III-22 

The  normalized  spectrum  for  heave  is  thus 


HN 


,  „  %  HhUfB.)  . 

“-Srcr 

Oh 


111-23 


Expanding  the  numerator  and  denominator  and  cancelling 

the  common  terms  S  (to)  and  |r_|*  this  leaves 

n  » 

HN  («,B,)  -  IXI-24 

h  I3h(») 

where 


nh(“'B*’  *  s  ¥•>  i?j* 


4  Cfo)  *•+*/* 

y5=  ~Z -  /  (1+cos (T  cos  B)  )  * 

/2*  B ,  -w/2 


(l+cos(T  sin  B)}  exp( *)  dB 

2o* 


13  (oj) 
n 


■W10 


sn<«>  |r#|* 


-  V/2 

(-)  /  (1+cos (T  cos  B) }{ 1+cos (T  sin  B))  dB 

*  -w/2 
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Similarly  for  roll 


Hr(u,$#)  II  (<u,0o) 
HNr(«,8.)  -  Hor(<u)  -  I3r(w) 


where 


II  («,8,> 
r 


Hr(w,8») 

s„U>  |r,|» 


d1  CM 
/lit  a 


8i+ir/2 

/  {1+cos (T  cos  8)) 
8,-ir/2 


X 


{l-cos(T  sin  0)}  exp{~^~rt~ 


13  (oi) 
r 


HQr 

sn(«) 


T 


d*  T/2 

-  /  { 1+cos (T  cos  8) > (1-cos (T  sin  0) >  d0 

11  -ir/2 


d8 
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IV.  Numerical  Prediction  of  Exciting  Forces 

A.  Discussion. 

Because  of  the  complicated  forms  of  the  integrands  of 
Equations  III-17  and  III-19  no  closed  form  solutions  to 
the  integrals  were  found  to  exist.  Numerical  methods  were 
employed  to  determine  the  variation  of  resulting  forces 
caused  by  the  variation  in  input  parameters.  Appendices 
I  and  XX  provide  program  listings  and  an  explanation  for 
their  use. 

B.  Assumed  values. 

In  order  to  use  numerical  methods  of  calculation 
values  must  be  assumed  for  the  relevant  parameters.  In 
this  thesis  the  initial  values  assumed  come  from  two 
Amoco  papers  published  in  the  "Journal  of  Petroleum 
Technology."  (See  Reference  1.)  These  papers  discuss 
a  proposed  vertically  moored  platform  (VMP)  which  is 
identical  to  the  concept  of  a  TLP.  Amoco 1 s  VMP  is  proposed 
to  have  leg  spacings  of  160  feet  and  is  estimated  to  have 
heave,  pitch,  and  roll  natural  periods  of  approximately  two 
seconds.  Because  of  the  popularity  of  the  cosine  squared 
wave  spectrum  a  comparable  standard  deviation  of  0.65 
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radians  was  chosen 


The  parameter  of  interest  in  all  of  the  following 
analyses  is  HU,  the  ratio  of  the  magnitude  of  the  transfer 
function  for  Gaussian  wave  spreading  to  the  magnitude  of 
the  transfer  function  for  uniform  wave  spreading. 

When  the  Amoco  TLP  was  analyzed  in  normally  distributed 
waves  with  a  0.65  standard  deviation  the  results  were  some¬ 
what  surprising.  The  variation  in  HN  from  unity  in  heave 
was  less  than  one  percent  for  all  values  of  wave  spectrum 
incidence/  B# ,  (see  Table  IV- 1) .  A  much  greater  variation 
of  HU  for  roll  was  expected  but  HU  still  remained  within 
ten  percent  of  unity  (see  Table  IV-2) .  Clearly  for  this 
case  an  initial  assumption  of  uniform  wave  spreading  and 
the  resulting  simplifications  in  calculation  of  the  force 
spectrum  will  provide  very  useful  first  pass  information 
to  the  designer. 

However  the  Amoco  proposal  is  by  no  means  the  only 
viable  configuration  for  a  TLP.  One  must  consider  the 
variation  in  HN  as  other  parameters  vary. 

C.  Variation  in  leg  spacing. 

The  most  obvious  variation  in  a  platform  is  the  spacing 
between  the  legs.  With  the  160  feet  of  the  Amoco  platform 
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as  a  reference,  leg  spacings  ranging  from  100  to  250  feet 
were  considered.  Analysis  results  for  heave  are  presented 
in  Table  IV- 3  and  for  roll  (and  hence  pitch)  in  Table  IV-4. 

For  all  of  these  configurations  HN  for  heave  is  still 
within  one  percent  of  unity.  Roll  suffers  somewhat  but 
the  calculated  HN  stays  well  within  20%  of  one.  For  initial 
design  work  this  accuracy  is  sufficient  and  the  assumption 
of  uniform  spreading  would  be  justifiable. 

D.  Variation  in  period. 

Because  the  wave  spectrum  is  continuous  in  frequency 
waves  of  all  periods  will  be  incident  upon  the  structure. 
However,  the  period  of  paramount  concern  is  the  natural 
period  of  the  structure.  This  varies  from  structure  to 
structure  and  is  somewhat  controllable.  Again  taking  the 
Amoco  value,  two  seconds,  as  a  reference,  periods  from  one 
to  three  seconds  were  considered.  The  leg  spacing  was 
assumed  to  be  160  feet  and  the  standard  deviation  0.65 
radians.  Analysis  results  for  heave  are  tabulated  in 
Table  IV-5.'  The  results  for  roll  however  vary  enough  to 
be  visible  on  a  polar  plot  of  HN  vs.  £• ,  the  center  of 
the  wave  spectrum,  and  are  shown  in  Figure  XV-6.  When 
possible  a  plot  has  been  preferred  to  a  table  because  of 


the  greater  amount  of  information  compactly  presented. 

For  heave  the  variation  in  HN  remains  at  about  one  per¬ 
cent  reaching  a  high  of  two  percent  for  a  three  second  per¬ 
iod.  For  roll  HN  shows  marked  fluctuation  but  again  stays 
within  twenty  percent  of  unity.  While  the  natural  period 
has  significant  implications  with  respect  to  the  periods  of 
significant  wave  energy  it  does  not  appear  to  significantly 
alter  the  modal  force  transfer  functions  as  long  as  some 
directional  spreading  is  present. 

E.  Variation  in  standard  diviation 

The  final  parameter  to  vary  is  the  standard  deviation. 
The  standard  deviation  for  uniform  wave  spreading  is  infinite 
while  that  for  unidirectional  waves  is  zero.  However,  as 
shown  in  Section  III-B  standard  deviations  above  two  radians 
approximate  the  condition  of  uniform  wave  spreading.  There¬ 
fore  almost  all  of  the  effect  of  varying  standard  deviation 
occurs  between  zero  and  two.  One  would  thus  expect  a  signi¬ 
ficant  change  in  HN  as  standard  deviation  changes,  particu¬ 
larly  for  small  values. 

This  expectation  is  correct.  Standard  deviation  was 
varied  from  0.125  to  2  radians  for  both  heave  and  roll 
(D  ■  160  feet.  Per  »  2.0  seconds).  The  results  are  plotted 
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in  Figures  IV- 7  and  IV-8  for  heave  and  roll  respectively. 

Consider  heave  (Figure  IV-7) .  For  standard  deviations 
of  1.0  or  2.0  radians  HN  stays  within  one  percent  of  unity. 
When  a  ■  0.5  the  variation  becomes  five  percent.  For 
o  ■  0.25  it  exceeds  twenty  percent  and  when  the  standard 
deviation  is  0.125  radians  even  the  pattern  of  the  curve 
changes.  For  roll  (Figure  IV-8)  the  results  are  similar, 
although  somewhat  more  pronounced. 

However,  if  one  can  allow  a  15%  error  in  an  initial  de¬ 
sign  approximation  then  the  uniform  wave  spreading  assumption 
is  valid  for  wave  directional  spectra  with  standard  deviations 
greater  than  approximately  0.5  radian.  Note  that  cosine 
squared  spreading  is  included  in  this  category. 

F.  Summary 

Numerous  other  cases  were  considered,  simultaneously 
varying  leg  spacing,  period,  and  standard  deviation.  In  all 
cases  where  o  >  0.5  the  resulting  normalized  amplitude- 
platform  force  transfer  function  had  a  magnitude  within 
20%  of  unity.  Thus  for  waves  with  any  reasonable  amount  of 
spreading  sufficient  accuracy  for  initial  design  calculations 
can  be  obtained  by  using  the  mean  square  value  of  r(a>,B). 

This  result  allows  for  a  simpler  initial  calculation  of 
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platform  response  to  random  waves  in  the  manner  of 
Vandiver.  For  situations  when  a  <  0.5  or  when  more 
accuracy  is  required  Vandiver's  method  can  still  be 
employed  using  the  computer  programs  of  this  thesis  to 
determine  HN  (Vandiver's  Cl) . 


SQUARE  TENSION  LEG  PLATFORM,  REAVE  EXCITATION 


WITH  GAUSSIAN  WAVE  DISTRIBUTION 


LEG  SPACING  -  160  feet 

PERIOD  ■  2.0  seconds 

STANDARD  DEVIATION  -  0.65  radians 


Correction 

factor  due  to  loss  of 

tails  -  1.016 

Factor  due 

to  uniform  spreading 

(13)  -  3.94 

a. 

11 

HNj^-11/13 

0 

3.90 

0.990 

5 

3.91 

0.992 

10 

3.92 

0.994 

15 

3.93 

0.997 

20 

3.94 

1.000 

25 

3.95 

1.002 

30 

3.96 

1.004 

35 

3.96 

1.006 

40 

3.96 

1.006 

45 

3.97 

1.007 

Table  IV-1  Normalised  Platform  Heave  Exciting  Force 
(D«160  feet,  Per-2.0  sec,  0*0.65  rad.) 
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SQUARE  TENSION  LEG  PLATFORM,  ROLL  EXCITATION 


WITH  GAUSSIAN  WAVE  DISTRIBUTION 


LEG  SPACING  -  160  £eet 

PERIOD  *  2.0  seconds 

STANDARD  DEVIATION  *0.65  radians 

Correction  factor  due  to  loss  of  tails  ■  1.016 
Factor  due  to  uniform  spreading  (13)  *  0.235  E+05 


8. 

I  1 

HNr-H/I3 

0 

0.216  E+95 

0.919 

5 

0.216  E+05 

0.920 

10 

0.217  E+05 

0.922 

15 

0.218  E+05 

0.925 

20 

0.219  E+05 

0.933 

25 

0.222  E+05 

0.943 

30 

0.224  E+05 

0.954 

35 

0.227  E+05 

0.967 

40 

0.231  E+05 

0.981 

45 

0.234  E+05 

0.996 

50 

0.238  E+05 

1.011 

55 

0.241  E+05 

1.026 

60 

0.244  E+05 

1.040 

65 

0.248  E+05 

1.053 

70 

0.251  E+05 

1.068 

75 

0.254  E+05 

1.079 

80 

0.255  E+05 

1.086 

85 

0.256  E+05 

1.090 

90 

0.257  E+05 

1.092 

Table  IV-2  Normalized  Platform  Roll  Exciting  Moment 
(D*160  feet,  Per»2.0  sec,  9*0.65  rad.) 


SQUABS  TENSION  LEG  PLATFORM,  REAVE  EXCITATION 
WITH  GAUSSIAN  WAVE  DISTRIBUTION 


PERIOD  •  2.0  seconds 

STANDARD  DEVIATION  -  0.65  radians 

CORRECTION  FACTOR  DUE  TO  LOSS  OF  TAILS  -  1.016 


LEG  SPACING 

HNh 

HNh 

(feet) 

(min.) 

(max. ) 

100 

0.997 

1.003 

150 

0.998 

1.003 

200 

0.994 

1.006 

250 

0.995 

1.004 

Table  IV- 3  Normalised  Platform  Heave  Exciting  Force 
(Various  leg  spacings,  Per-2.0  sec,  O-0.65  rad.) 
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SQUARE  TENSION  LEG  PLATFORM,  ROLL  EXCITATION 
WITH  GAUSSIAN  WAVE  DISTRIBUTION 


PERIOD  *  2.0  seconds 

STANDARD  DEVIATION  -  0.65  radians 

CORRECTION  FACTOR  DUE  TO  LOSS  OP  TAILS  -  1.016 


LEG  SPACING 

HNr 

HNr 

(feet) 

(min.) 

(max.) 

100 

0.978 

1.027 

150 

0.933 

1.067 

200 

0.882 

1.114 

250 

0.846 

1.160 

Table  IV-4  Normalised  Platform  Roll  Exciting  Moment 
(Various  leg  spacings.  Per* 2.0  sec,  o*0.65  rad.) 
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SQUARE  TENSION  LEG  PLATFORM,  HEAVE  EXCITATION 
WITH  GAUSSIAN  WAVE  DISTRIBUTION 


LEG  SPACING  - 
STANDARD  DEVIATION  - 

CORRECTION  FACTOR  DUE 


PERIOD 

HNh 

(seconds) 

(min.) 

1.0 

0.997 

1.5 

0.996 

2.0 

0.990 

2.5 

0.994 

3.0 

0.980 

160  feet 
0.65  radians 

TO  LOSS  OF  TAILS  -  1.016 

(max.) 

1.001 

1.005 

1.007 

1.004 

1.024 


Table  IV- 5  Normalized  Platform  Heave  Exciting  Force 
(D«160,  Various  periods,  a-0.65  rad.) 
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Figure  IV-6  Normalized  Platform  Roll  Exciting  Moment 
(1>160,  various  periods,  c-0.65  rad.) 
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Lgure  IV- 7  Normalized  Platform  Reave  Exciting  Force 
(D-160,  Per  -  2.0  aec,  varioua  a) 
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V.  Suggestions  for  Experimental  Verification 

A.  Single  leg. 

Because  the  results  of  the  computer  program  provide  a 
transfer  function  between  forces  on  a  single  leg  to  platform 
forces  and  moments,  the  single  leg  forces  must  be  obtained. 
This  can  be  done  theoretically  or  experimentally.  This 
report  will  deal  with  experimental  determination. 

The  single  leg  is  assumed  to  be  of  circular  cross  sec¬ 
tion.  Thus  there  is  no  angular  dependence  and  testing  can 
be  performed  in  a  ship  model  towing  tank  with  a  wave  gener¬ 
ator.  In  the  proposed  test  the  single  leg  model  will  be  held 
rigid  with  load  cells  measuring  the  forces  as  a  result  of 
incident  waves.  These  forces  will  be  normalized  to  forces 
due  to  scaled  unit  amplitude  waves  using  a  linear  assumption 
for  the  force  to  wave  amplitude  relation.  By  performing 
this  test  with  a  variety  of  wave  frequencies  a  plot  of  rs(a») 
can  be  obtained.  From  this  the  computer  will  predict 
platform  heave  force  and  roll  moment  transfer  functions. 

B.  Square  TLP. 

To  completely  test  the  results  of  this  study  would  re¬ 
quire  the  ability  to  generate  a  directional  wave  spectrum. 
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This  ability  is  not  available.  However,  the  total  transfer 
function,  rt(w,3) ,  can  be  experimentally  verified  for  heave, 
pitch,  and  roll. 

Once  the  results  for  a  single  leg  are  obtained  a 
similar  experiment  with  the  total  platform  can  be  performed. 
The  model  should  be  held  rigid  with  load  cells  measuring 
forces  in  all  six  degrees  of  freedom.  The  direction  of 
wave  incidence  can  be  varied  by  turning  the  model. 

This  experiment  will  generate  a  family  of  curves  of 
rt(o»,3)  ,  one  curve  for  heave  and  one  for  roll  for  each  angle 
of  incidence  3.  These  curves  can  then  be  compared  to  those 
predicted  using  the  measured  forces  for  rg(o>)  and  the  heave 
and  roll  algorithms  developed  in  this  paper.  Equations  III-7 
and  III- 10. 
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VI .  Conclusions 


A.  Assumptions  and  approximations. 

1.  Linearity.  The  forces  computed  in  this  paper  are 
those  exerted  by  waves  on  a  motionless  structure.  The 
forces  are  assumed  to  be  linearly  related  to  wave  amplitude. 
Furthermore  the  forces  exerted  on  a  motionless  structure 
can  only  be  useful  in  predicting  dynamic  response  if  the 
equations  of  motion  are  linearized.  For  a  TLP  the  equations 
of  motion  will  be  linear  for  low  to  moderate  sea  states. 

This  is  due  to  the  large  dimensions  of  the  major  members  and 
to  the  relatively  small  response  amplitudes.  Damping  terms 
for  the  TLP  may  also  be  linearized  because  the  damping 
ratios  are  expected  to  be  small  in  heave,  pitch,  and  roll. 

2.  Effects  of  cross  members.  The  model  of  a  TLP  used 
to  derive  the  transfer  functions  had  four  legs  acted  upon  by 
wave  forces,  and  nothing  else.  Real  TLP's  must  have  cross 
members  between  the  floatation  legs  to  provide  structural 
support.  These  would  also  be  subject  to  wave  forces. 

The  diameter  of  a  cross  member  will  be  much  smaller  than 
that  of  the  leg.  In  the  Amoco  platform  referred  to  in  Refer* 
ence  1  the  cross  members  have  an  outer  diameter  of  5.5  feet 
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while  the  floatation  leg's  diameter  is  30  feet  near  the  sur¬ 
face  where  the  large  wave  forces  are.  This  is  a  factor  of 
about  5h  difference.  The  displaced  volume  is  proportional 
to  the  cross  sectional  areas  which  are  almost  a  factor  of  30 
different.  In  addition  the  support  members  are  mostly  near 
the  bottom  of  the  legs,  away  from  the  largest  wave  forces. 
The  parts  that  are  shallow  are  close  to  the  legs  and  would 
react  in  phase  with  them.  Thus  the  effect  of  the  waves  on 
the  cross  members  can  be  assumed  to  be  negligible  with 
respect  to  the  effect  on  the  legs. 

3.  Wave  disturbances  of  the  legs.  Another  aspect  not 
considered  in  this  model  is  the  effect  of  one  leg  upon  the 
waves  incident  upon  the  others.  In  the  Amoco  example  the 
outer  edges  of  the  legs  at  the  water  surface  are  130  feet 
apart  which  is  only  about  four  times  their  30  foot  diameters. 
There  could  be  some  shielding  or  reflecting  effects  present, 
causing  the  results  of  this  study  to  be  in  error.  As  the 
ratio  of  leg  spacing  to  leg  diameter  increases  these  errors 
will  decrease. 

The  primary  purpose  of  this  thesis  was  to  evaluate  the 
influence  of  wave  spreading  on  the  heave,  pitch,  and  roll 
force  spectra.  The  most  important  geometric  effect  is  the 
relative  phases  of  forces  exerted  on  the  legs.  In  cases 
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that  the  wave  spreading  is  shown  to  be  sufficient  to  average 
out  the  relative  phase  effects,  then  it  can  certainly  be 
argued  that  interaction  effects  between  legs  would  also 
be  suppressed. 

B.  Utility  of  results. 

1.  Initial  design  approximation.  For  wave  directional 
spreading  with  a  standard  deviation  of  greater  than  0.5 
radians  the  heave  force  or  roll  or  pitch  moment  imposed 
upon  a  TLP  can  be  approximated  within  20%  to  be  that  imposed 
by  uniformly  spread  seas.  By  incorporating  this  into  the 

4 

method  used  by  Vandiver  an  initial  estimate  of  the  response 
of  a  TLP  can  be  obtained  which  would  be  applicable  in 
preliminary  design. 

2.  Mere  accurate  approximation.  If  the  standard 
deviation  of  the  wave  spectrum  is  less  than  0.5  radian  or 
if  an  answer  with  more  than  20%  accuracy  is  required  then 
the  computer  program  outlined  in  Appendices  I  and  II  can 
be  used  to  determine  HH  (Vandiver's  Cl).  This  can  then 

be  incorporated  into  the  calculations  of  the  modal  response 
of  the  structure. 

3.  Variation  to  fit  other  structures.  Results  for 


other  square  TLP's  can  easily  be  obtained  simply  by  changing 


the  values  assigned  for  the  leg  spacing  and  the  period  of 
interest.  If  the  structure  of  interest  is  not  square  the 
same  concept  can  be  applied  by  changing  the  program  defini¬ 
tions  of  the  transfer  functions.  In  this  manner  results 
can  be  obtained  for  any  shape  of  TLP  (square,  triangular, 
rectangular,  etc.)  with  any  leg  spacing  and  any  natural 
period. 

4.  Wave  directional  spectrum  variation.  The  greatest 
advantage  of  modeling  the  wave  directional  spectrum  as  a 
Gaussian  process  is  the  continuous  variation  in  spreading 
it  provides.  The  standard  deviation  of  the  wave  spectrum 
at  a  proposed  location  for  a  TLP  could  be  measured  under 
various  conditions,  noting  that  the  periods  of  the  signi¬ 
ficant  wave  energy  would  also  change.  This  information, 
standard  deviation  and  period,  could  then  be  input  to  this 
program  to  estimate  platform  forces  and  moments  under  those 
conditions.  The  final  information  of  forces,  moments,  and 
responses  would  then  be  applied  to  sizing  the  tension 
members  for  adequate  tensile  and  fatigue  strength  and  to 
ensure  the  platform  motions  do  not  cause  equipment  malfunc¬ 
tions  or  personnel  discomfort. 
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NOMENCLATURE 


a  distance  from  center  of  TLP  to  center  of  any  leg 

C  constant  to  normalize  the  area  of  the  wave 

directional  spectrum 

d  leg  spacing  of  TLP 

g  acceleration  of  gravity 

h  as  a  subscript,  denotes  heave 

Hx  frequency  spectrum  of  the  total  force  or  moment 
on  the  TLP  due  to  wave  spreading 

Hqx  frequency  spectrum  of  the  total  force  or  moment 

on  the  TLP  due  to  uniformly  distributed  waves 

HNX  normalized  frequency  spectrum  of  the  total  force 
or  moment  on  the  TLP 

Tx  integral  portion  of  Hx  dependent  on  wave  spreading 

Ilx  numerator  of  HNX,  dependent  on  platform  geometry 
and  wave  spreading 

I3X  denominator  of  HNX,  dependent  only  on  platform 
geometry 

k  wavenumber 

1^  distance  from  center  of  TLP  to  projection  of  leg  i 
on  the  direction  of  wave  propagation 

p  as  a  subscript,  denotes  pitch 

r  as  a  subscript,  denotes  roll 

R  ratio  of  magnitude  squared  of  platform  r  to 
magnitude  squared  of  single  leg  T 

s  as  a  subscript,  denotes  single  leg 

Sn  wave  spectrum 
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t  as  a  subscript,  denotes  total  platform 

T  non-dimensional  quantity  determined  by  leg  spacing 
and  wave  frequency 

6  angle  of  wave  incidence  measured  from  roll  axis 

St  angle  of  incidence  of  center  of  wave  spectrum 

r  transfer  function  from  incident  wave  amplitude  to 

x  force  (or  moment)  exerted  upon  the  structure 

5„  Dirac  delta  function 

inn 

X  wavelength  of  incident  waves 

a  standard  deviation  of  wave  directional  spectrum 

E  summation 

phase  shift  of  the  wave  incident  on  leg  i  as  compared 
to  the  wave  at  the  center  of  the  TLP 

*  normal  distribution  function 

o)  radian  frequency  of  incident  waves 
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Appendix  I .  Computer  Program  Listing 

A.  Discussion. 

Because  of  the  complicated  forms  of  the  integrands  of 
Equations  III-17  and  111-19  no  closed  form  solution  to  the 
integral  was  found  to  exist.  Yet  the  trigonometric  func¬ 
tions  with  trigonometric  arguments  caused  such  rapid  oscil¬ 
lations  of  the  integrand  that  most  numerical  methods  failed 
to  converge.  The  method  found  to  be  the  most  useful  was 
Gauss-Legendre  quadrature  (see  Appendix  III) .  A  quadrature 
of  fifth  order  was  found  to  give  sufficient  numerical 
accuracy  and  is  contained  in  the  subroutine  GAULEG. 

Because  the  wave  directional  spectrum  was  only  non¬ 
zero  between  plus  and  minus  ir/2  radians  a  correction  factor 
was  necessary  to  ensure  its  area  was  unity.  This  correction 
was  accomplished  by  an  eighth  order  polynomial  approximation 
to  the  area  under  the  normal  curve.  This  approximation  was 
developed  by  the  author  and  is  contained  in  the  subroutine 
ERF. 

Because  of  peculiarities  in  the  computer  system 
available  the  platform  transfer  functions  are  entered  as 
function  subprograms  HEAV  and  FKEAVE  for  heave,  and  ROLL 
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and  FRQLL  for  roll,  the  former  of  each  pair  being  the 
transfer  function  with  uniform  wave  spreading. 

The  program  is  written  in  FORTRAN-IV  language  and  has 
numerous  comment  statements  to  aid  in  its  understanding. 

B.  Program  listing. 

A  listing  of  the  programs  for  heave  and  roll,  with 
associated  subprograms,  follows. 


[ 

! 
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TESTIS  IE3  PLATFOM 
«?«»! 

ESTEPIM  PSAT,Ff»EAf 

SET  ST**OUi>  DFTIATIO*  (3D), 155  SPACHC  (9), AID  ?E»rOD  (PER) 

0«153. 

S9-0.55 

PF9«2. 

CUC0LA Tl  VAIIAUCE  *19  T 
293  m»59**2 

333  T»1.226:35*0/*ER«*2 

V»I*E  ( 5,5) 

5  s'S’-AT  (*1*,T1'*, ‘SQUARE  TF1SI31  tS3  PL*  TF0K,*,B6P9S  EXCmTX39‘> 
WPITC  (5,1?) 

19  ?9*1».r  C  ‘.TII.UACSSTAI  WA7E  DISTRIBUTim*//) 

W»m  (5,13)  3 

15  FOP*»T  (‘  *  ,*15  ,*  LES  5P»Cri5»*,Ttt3,F1.3,  *  FEET* ) 

SPITE  (5,23)  PF» 

23  *,3°,'1T  C  •,ri«,‘PEPm»‘,?0?,*1.1,‘  SEC3P9F* ) 

WHIT?  (6,25)  59 

25  *9»«H  C  ‘,ri«,‘ST"!9A»D  0EU*rr3l«‘,T?1,F'.2,‘  RADIARSV/) 

cappscrroi  fast?*  f9»  loss  ar  mt?  of  saussiai 
ar9T1«1.*!’?7n5/S9 
cut  ERF  (ai9M,P1I) 

C»1./(2.*P“I-1.) 

»*?**  (5,13)  : 

33  FCP’'**  (•  ‘.riS.'C*'*  RECTI  OS  FACT3P  30E  T3  13SS  3F  T*ItS«‘  ,F6.3) 
FACTOR  905  T9  UUF98F  EPSFSOrP'J 

Et»-1. 579796 
T3»  1.570796 

cut  3AUVC  (U,X3,HUT,R,IP,»-*,F*F,T> 

Rt3*R 

WRITE  (5,15)  111 

35  F9P-«T  C  •,T1«,,F»Cr9®  ODE  -3  nilFOPH  ?®RE*DT1S  (I3)»‘,E11.*> 

mv  (5. »?)  r* 

a?  F3?.'*AT  C  *  »T2A»  *  (EHP3*  CaRrillO***  ,13,  *)  V) 

(5,05) 

»5  F9**AT  (•  •,*  B?r*,,r23,,ri,,Ttt5,,E9»I1/ri,,T5?,,SRSORV/) 

VA’T  9IFECTI31  3F  CESTE1  O'  HAYS  STECTP9* 

103.  23»830F0‘1.7»5329E-2 

U1»Et*5? 
r»1«YU*5C 

CALL  2AUIE3  (U%<01,R9E>»,T,IEP.,10,T*F,T) 

8T1«T*C/S9 

H1«iri/(*I3 

WRITE  (6,53)  899S'5,9Ti,*9,IEP 
59  n***T  (•  •,F5.1,"17,S11.0,T02,F10,»,T52,I1) 

■5P“»«B39P3*5. 

IF  (P30E3  .IF.  05.)  90  TO  190 
5T3P 
£10 


THIS  PAGE  IS  BEST  QUA1ITT  PKACTICABLS 

Cvi  )  L  x  'O  1,0'C 
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U  U  o  u  u  u  u  u  o  u 


I 


TENSION  LEG  FLATFORP 
NOLL 

EXTERNAL  ROLL .FROLL 

SET  STANQARO  DEVIATION  (SO).LfG  SPACING  (0  )  »  A  NC  PERIOD  (PER) 

0*160. 

SO*G .69 
PER*?. 

CALCULATE  VARIANCE  ANO  T 
2®C  VAP*SD**2  • 

300  T*1.22603«*0/PCR*«2 

WRITE  (6.S ) 

3  PORPAT  (*1*»T10»*S6U*RE  TENSION  LEG  PL'TFORN.  ROLL  EXCITATION*) 
WRITE  (6.10) 

10  FORPAT  (*  • »TIS. 'GAUSSIAN  WAVE  DISTRIBUTION*//) 

WRITE  (6. IS)  0 

IS  FORMAT  »*  '.T13.  *LEO  SPACING* • *TA O.FS. 0 ••  FEET*) 

WRITE  (6.20)  PER 

2C  FORMAT  (•  *.T15.*PERIOr**.T*O.FS.l**  SECONDS*) 

WRITE  (6.25)  SO 

23  FORPAT  (»  •«T13.*STAX0AR0  DEVIATION*' *  T39.F6.2.'  RADIANS'//) 

CORRECTION  FACTOR  FOR  LOSS  OF  TAILS  OF  GAJSSIAN 
WinTH*1.5707?6/SO 
CALL  ERF  (WIDTH, PHI) 

C*l./(2. .PHI-1.) 

WRITE  (6.30)  C 

30  FORPAT  (•  ',T15.«C0RRECTIPN  FACTOR  DUE  TO  LOSS  OF  TAILS** »FG.3 ) 
FACTOR  DUE  TO  UNIFORM  SPREADING 

XL*-i. 370796 
XU*  1.370796 
57*0. 

call  gaulen  ( xl« xu.ro ll • r. 1 r« a0. VAR. T) 

RI3*R'0*»2 
WRITE  (6.33)  RI3 

33  FORPAT  (•  *.T13. 'FACTOR  DUE  TO  UNIFORP  SPREADING  ( 13) *• .Ell .4 ) 
WR1TC  C6.*C)  IR 

♦0  FORPAT  <•  • «T20 .* (CRPOR  CONDI TION* •» IS /) 

WRITE  (6«*5) 

♦9  FORPAT  (•  BETA  *» 72 1.*Il*tT*5**MN*II/I3*»T62»* ERROR?*/) 

VARY  DIRECTION  OF  CENTER  OF  WAVE  5PCCTRUP 
BOOEGSO. 

100  R0*"0CEP*1.7A3S2RE-? 

Xll*XL«93 

XU1*XU*D0 

call  oauleo  (xll.xul . Troll *y.I or. 30. VAR.T) 

A!1*Y.C**2*C/SD 

HRSRI1/RIS 

WRITE  (6*30)  6’IOCG.R  I  l.PN  «TCR 
30  FORPAT  (•  **FS. 1.717. CU.*.TA2.F1Q.«.T6?.I3) 
tOOEG*BQDEG«f. 

IF  (ROOEG  .LE.  90.)  GO  TO  100 

STGP 

CNO 


JJJIS  PACT 13  ' 

Wj*  hv 


jAiiin 

,  x  J  u>C 
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±  *■ 


so??oorivr  (it,nj,pc?,T,T?<>,o,),Y?»,T) 
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FUNCTION  FHEAV(3*5D*V'AR  »T> 

FHEAV  =  l.f‘»*>76“*<l«*C  n-t*  T*COf  )  >  5  »  <  l.-C  E  *  T  *$  IM(  •?  ‘  >  )  • 
1  EX?C-<;i-no>**2/<2.*VAr)) 

RETURN 

END 


function  HE  AV  {b*  ?  0  •  V 

HEAV  =  1.27324*»1«-»C0S<T*CCS(B))  )*(  1  .  ♦ 1 E  <  T  *  S  I  « *  3  *  )  > 

RETURN 

END 
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A.  Discussion 


This  program  was  intentionally  designed  to  remain  as 
much  as  possible  in  discrete  sections,  each  performing  a 
specific  purpose.  In  this  manner  the  program  can  be  more 
easily  understood  by  others  and  converted  to  use  for  many 
different  sizes  and  types  of  structures. 

An  error  flag  has  been  added  to  the  subprogram  GAULEG. 

If  GAULEG  is  unable  in  a  reasonable  time  to  get  consecutive 
iterations  to  agree  within  one  percent  it  flags  this  fact 
to  the  user  by  indicating  an  error  condition  of  one.  If  the 
error  condition  is  zero  then  GAULAG  has  gotten  the  same  an¬ 
swer  (within  one  percent)  for  two  iterations  with  the 
second  iteration  having  twice  as  many  increments  as  the  first. 

B.  Numerical  parameters. 

To  make  it  easier  to  understand  the  program  flow  the 
leg  spacing  (D) ,  standard  deviation  (SD) ,  and  period  (PER) 
are  simply  assigned  values  in  FORTRAN  statements.  They  could 
just  as  easily  be  read  from  data  cards,  input  from  an  inter¬ 
active  terminal,  or  varied  in  do-loops,  whichever  is  most 
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convenient  for  the  user  on  the  available  system. 


To  adopt  the  program  to  a  square  TLP  of  interest  one 
must  simply  replace  the  assignment  statements  for  leg  spacing 
and  for  the  period  of  interest.  If  many  periods  are  import¬ 
ant  then  a  do- loop  for  period  should  be  used  returning  to 
statement  300  to  ensure  the  value  of  T  is  updated.  If  a 
do-loop  is  used  to  vary  standard  deviation  it  must  return 
to  statement  200  to  revise  the  value  of  the  varience. 

The  center  of  the  wave  spectrum  ( 3  o )  was  varied  in  an 
implied  do-loop  with  statement  100  because  one  is  normally 
interested  in  the  response  to  waves  coming  from  all  directions. 
Because  of  the  symmetry  properties  of  a  square  TLP  the  angle 
must  only  be  varied  over  45  degrees  for  heave  and  90  degrees 
for  roll. 

C.  Conversion  to  other  types  of  platforms. 

The  information  concerning  the  platform  geometry  is 
contained  in  the  platform  transfer  functions.  In  the  computer 
program  the  definition  of  the  functions  in  the  function 
subprograms  (ROLL,  FROLL ,  HEAV,  and  PHEAV)  entirely  defines 
the  structure.  Note  that  a  factor  of  d*  was  left  in  the  main 
program  rather  than  ROLL  and  FROLL  since  this  factor  will  be 
present  in  moment  transfer  functions  for  any  type  of  structure. 
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Thus  to  change  this  program  to  accomodate  any  shape  TLP 
one  must  simply  perform  derivations  similar  to  those  in  sec¬ 
tion  III-A  and  use  the  results  to  replace  the  function  defin¬ 
itions  in  the  function  subprograms. 

D.  Conversion  to  other  wave  spectra. 

This  program  was  explicitly  written  for  a  normally  dis¬ 
tributed  wave  directional  spectrum  because  it  is  so  easily 
varied.  Simply  by  varying  the  standard  deviation  one  can 
obtain  as  directional  or  uniform  a  spreading  as  desired. 
However,  if  another  spectrum  is  preferred  this  program  can 
be  adapted. 

The  most  obvious  place  to  change  is  in  the  definition 
of  FROLL  and  FHEAV.  There  the  new  spectrum  would  be  substi¬ 
tuted  for  the  exponential  term.  In  the  main  program  the 
correction  factor  C/SD  is  applied  after  calling  GAULEG  using 
FROLL  or  FHEAV.  This  would  be  changed  or  eliminated.  Since 
C,  the  factor  for  loss  of  tails,  is  no  longer  necessary  the 
entire  subroutine  ERF  could  then  be  eliminated.  New  correc¬ 
tion  factors  for  the  new  spectrum  would  then  be  introduced. 

It  is  also  possible  to  normalize  the  transfer  function 
using  something  other  than  uniform  spreading.  For  this  ROLL 
and  HEAV  only  must  be  revised.  If  correction  factors  are 


65 


necessary  they  are  best  introduced  in  the  main  program  after 
calling  GAULEG  for  the  integration. 
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Appendix  III.  Gauss-Leqendre  Quadrature 


A.  Legendre  Polynomials 

1.  Definition.  Define  a  set  of  polynomials  Pn(x)  of 
degree  n,  where  n  is  a  non-negative  integer,  such  that 

Pn^x^  "  cn  **mn 

and  P_(l)  ■  1. 

n 


The  first  is  an  orthogonality  condition  and  the  second 
is  necessary  to  fix  cn.  Together  they  define  a  unique  set  of 
polynomials  known  as  the  Legendre  polynomials,  the  first  few 
of  which  are 


PQ(x)  »  1 
P^x)  »  x 
P2(x)  «  ij(3x*-l) 

P3 (x)  «  H(5x*-3x). 


In  general  they  can  be  found  by 


P«(j°  "  2n  n|  ^n  Cx*-l)n. 


Note  that  P„(x)  has  n  seroes  in  the  interval  from  minus  one 
n 

to  one  and  is  even  for  n  even  and  odd  for  n  odd. 
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.  Polynomial  approximation.  Any  function  piecewise 


continuous  in  the  interval  from  minus  one  to  one  can  be 
expanded  in  terms  of  an  infinite  series  of  Legendre  polynomials 
It  can  be  approximated  by  a  finite  series  of  Legendre  poly¬ 
nomials  . 

yU)  "J.6®  P®<X)  %J,C*  P»(x> 

B.  Integral  approximation 

A  simple  means  to  evaluate  the  integral  in  the  region 
minus  one  to  one  is  as  follows: 

:  i  • 

I  y(x)  dx  ■  /  I  c-  Pm(x)  dx 

-i  -i  •  " 

Multiply  both  sides  by  1  *  PQ(x): 

}  y(x)  dx  m  Z  c  }  P_(x)  P  (x)  dx 
^  I  |  nu»  |  in  o 

-  c0  /  PQ(X)  PQ(x)  dx 
y(x)  dx  *  2  cQ 

Thus  the  problem  has  now  become  one  of  finding  cq»  the 
coefficient  of  PQ(x) .  One  possible  method  will  be  shown 
using  the  example  where  n-2. 

Me  desire  to  find  coefficients  and  arguments  such  that 
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A 


-I  y<*>  dx  -  1  C,  -  A,  y(x#)  +  ylx^  . 

For  any  two  points.  x#  and  x,  ,  we  expect 

C#  P,  (x,)  ♦  Ct  P^x,)  ♦  C2  P^x,)  -  y(X|) 

Ct  P.U,)  +  Cx  PjfXj)  ♦  C2  PjCx^  -  y  (x^ ) 

Noting  that  P  {x>  ■  1  and  P  (x)  »  x  and  choosing  x.  and 

■  l  0 

Xj  to  be  the  zeroes  of  P{ (x)  this  becomes 

C,  +  Ct  x#  -  y(x#) 

C,  ♦  Cj  xt  ■  y(x2) 

Where  x(  and  xl  are  t  (1//3)  .  By  adding  these  two 
equations  and  noting  that  x(  «*  -xt  we  have 

2  C,  -  y  (-1//3)  +  y(l//3) 

which  is  our  desired  answer. 

To  get  a  larger  degree  of  accuracy  a  larger  n  is  required. 
This  method  for  determining  2  C,  remains  effective.  For  the 
case  of  n  ■  5  the  solution  is  (to  seven  significant  digits)3 

i 

{y(x)  dx  -  (0.568  888  9)  (y(0>) 

♦  (0.478  628  7)  {y(.538  469  3)+y(-.538  469  3)} 

♦  (0.236  926  9)  (y(.906  179  8)+y(-.906  179  8)}  . 
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However,  one  is  not  always  interested  in  the  interval 
from  minus  one  to  one.  A  transform  must  then  be  used  to 
transfer  the  region  of  interest  to  that  interval.  For  addi¬ 
tional  accuracy  without  the  need  for  a  larger  order  quadra¬ 
ture  the  total  interval  can  be  divided  into  subintervals. 
Bach  subinterval  can  be  transformed  seperately  and  the 
result  summed.  Thus  each  subinterval  will  be  approximated 
by  a  fifth  order  polynomial.  This  is  the  method  employed 
in  the  subroutine  GAOLEG. 


